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We describe an efficient scheme for evaluating higher order contributions to primordial cosmo- 
logical perturbations using the "in-in" formalism, which is the basis of modern calculations of 
. . . non-Gaussian and higher order contributions to the primordial spectrum. We show that diagrams 

0^ ' with two or more vertices require careful handling. We present an implementation of the operator 

. formalism in which these diagrams can be evaluated in a simple and transparent fashion. We illus- 

' trate our methodology by evaluating the correction to the primordial gravitational wave spectrum 

04 ' generated by scalar loops, a 2-vertex, 1-loop interaction. We then look at a generalized A''-point, 

j^jQ, 2-vertex diagram. 
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I. INTRODUCTION 



The pairing of the Arnowitt-Deser-Misner (ADM) formulation of general relativity [l| and the "in-in" approach 
to quantum field theory [3, 0iJE| is a remarkably powerful tool for analyzing cosmological perturbations^. First 
^ app lied to cosmology by Jordan Q and Calzetta and Hu [9|], Maldecena's treatment of primordial non-Gaussianities 

I ' established it as the preferred approach to computing both iV-point correlation terms [ll|, 113, [Hj HI, [H, d, [i3| 

and loop corrections [H, [l^, [13, [2ll |22|. [23l . [2^ . [25| to the inflationary perturbation spectrum. The perturbations 
(— I ■ generated by very simple models of inflation are accurately described by the long established lowest-order expressions 
[2^ 113, [111 . However, current [1^ and forthcoming experiments [sO] will put tight limits on the primordial bispectrum 
(and possibly the trispectrum) , and thus place tight constraints on inflationary models with non-trivial perturbation 
spectra, so these calculations are of substantial practical importance. 

The purpose of this paper is to present an efflcient and transparent scheme for implementing the operator formalism 
[l8| . Additionally, we warn that in some formulations, diagrams with more than one vertex^ can become pathological if 
, their external momenta or vertices are not distinguished from one another in some way. The problem can be viewed as 

■ one of initial conditions - "in-in" takes the initial states of some set of fields, |in), and calculates the expectation value 
of some set of operators with respect to these fields at a later time t, \^{t)). Specifying the initial conditions typically 
amounts to choosing the Bunch-Davies vacuum: the dynamical degrees of freedom behave like harmonic oscillators at 

' very small length scales, and each mode is assumed to be in its ground state. There is then an operational prescription 

■ that enforces this selection in field theoretic calculations. As we show below, specifying the initial conditions implicitly 
restricts the algebraic manipulations one may perform. To calculate an expectation value at some arbitrary time in 
the interaction picture we need to evolve the in-state forward in time. Operationally, the simplest way of incorporating 
this initial condition is to include some evolution in imaginary time, t — > t{l + ie). There are (at least) two simple 
ways of expressing this choice. One may define the integration in the time evolution operator to run along a contour 
in the complex plane, rather than along the real axis. Alternatively one may analytically continue the time variable 
itself to include a small imaginary component, leaving the integral on the real axis. Importantly, some algebraic 
manipulations which are permitted in one formulation are forbidden in the other. 

Interestingly, this problem cannot arise in conventional "in-out" calculations, or in "in-in" calculations for diagrams 
with a single vertex - this discussion is thus timely, since theoretical analyses of primordial perturbations have matured 
to the point where more complicated interactions are now routinely considered [3lll32|. We use the operator formalism 
introduced by Weinberg ^18j , as it provides an efficient and transparent scheme for performing in-in calculations, and 
the specific approach we develop here can significantly reduce the algebraic workload associated with a given diagram. 

This paper is organized as follows. We review the operator formalism in Section [TTl using the notation summarized 
in the Appendices of our previous paper [23| . We review two different prescriptions for injecting the Bunch-Davies 



^ The formalism is the subject of two recent comprehensive monologues 0,0], where the reader is referred to for further information. 
^ We are not restricting our attention to diagrams with a fixed number of external legs, so the order of a diagram here is simply the 

number of vertices, and not the number of loops — the vertex-count is equivalent to the number of copies of the interaction Hamiltonian 

being correlated. 
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initial conditions, showing how to obtain consistent results within the formalism of Weinberg [18|. In the two following 
Sections we present sample calculations. In mil] we calculate scalar loop corrections to the graviton power spectrum. 
As one would expect, the scalar corrections to graviton (tensor) perturbations are too small to have any practical 
impact, but the calculation is a useful example of our overall methodology. In mV\ we write down a two vertex 
contribution to the iV-pt correlation function with p internal lines. With N ^ 2, p ~ 2 this is a correction to the 
propagator with the same topology as the graviton loop correction of ^ 3111) and for 4, p = 1 we have the topology 
of the graviton [3l| and scalar (sj exchange contributions to the scalar 4-point (trispectrum). We show that these 
diagrams can become problematic in limits where the external momenta are not distinguishable. We demonstrate 
that explicitly injecting the initial conditions by deforming the time contours into the complex plane at the outset 
sidesteps any difhculties and substantially simplifies the algebra. We conclude in Section FVl 



II. SPECIFYING INITIAL CONDITIONS IN "IN-IN' 



In the in-in formalism, calculations of the expectation value, {W{t)), of a product of operators W{t) at time t, 
require that we evaluate 



{W{t)) =(^(Te-'^-^"-'^''^''y W{t) (Te-'^-=.«..«(*")'^*")^ , 



(1) 



where the fields on the right hand side are Heisenberg fields. The expectation value is taken with respect to the initial 
state, |in), which we assume to be the Bunch-Davies vacuum. The interaction Hamiltonian Hint is defined in the 
usual way, so that the total Hamiltonian H is the combination of Hi^t and the free- field Hamiltonian, Hq, 

H^Ho + Hn,f (2) 

Analyses of cosmological perturbations typically start with the Einstein- Hilbert action of general relativity together 
with the appropriate matter action. One then uses the ADM formulation [l[ to obtain an action containing only 
dynamical degrees of freedom [l0|,[13]- From the action one constructs the Hamiltonian by defining conjugate mo- 
menta, and separating out the quadratic from the higher order parts: iJo consists of terms that are quadratic in the 
perturbative degrees of freedom (and thus free) , while iJi„t consists of all third and higher order terms [l^ . The free 
Hamiltonian Hq drives the evolution of the operators, while TJjnt evolves the states. This separation is natural, since 
in a homogenous and isotropic background we can find the eigenstates of the free field Hamiltonian at past infinity. 

The interaction terms generally have derivative couplings even when the action contains only canonical kinetic 
terms. These derivative couplings are the end result of perturbatively expanding the action. In models with non- 
standard kinetic terms such as DBI inflation ^33i] or k-inflation 34] , derivative couplings are generically present at the 
outset. Consequently, we must proceed carefully when canonically quantizing the theory. If Lint is the portion of the 
action with terms third order and higher, the usual expression for the interaction Hamiltonian i?i„t = —Lmt acquires 
corrections at fourth order, as first shown in [2^. In slow roll inflation these corrections are proportional to powers 
of the inflationary slow-roll parameters, but are generically unsuppressed in DBI or k-inflation [s^, [s^. In the more 
general case, for instance if one wants to calculate correlations in the radiation dominated era, there will be extra 
interaction terms. 

To make contact with the familiar "in-out" formalism of QFT, insert complete sets of states labeled by a and (3 
into equation ([1]), 

{W{t))= j da J dP{0\ (Te"'^*o^'°'^*'^''*')V)(a|I^(t)|/3)(/3| (re"*^'o ^■"'^*"^''*") |0) . (3) 

The interpretation is clear, the "in-in" correlation is the product of vacuum transition amplitudes ("in-out") and 
a matrix element (a|M^(t)|/3), summed over all possible "out" states. The "in-in" formalism is simply standard 
QFT, rigged to compute correlation functions at fixed time, given initial conditions instead of asymptotic boundary 
conditions. Initial conditions in QFT are usually specified by finding the eigenstates of the free Hamiltonian Hq, 
and stipulating that the system begins in one (or some combination) of these eigenstates. If the system begins in 
the quantum mechanical vacuum, this amounts to putting our system in the vacuum state of Hq at the initial time. 
Operationally, the vacuum is selected by redefining the range of t to include a small imaginary component [3^ . 
t ~>t + i€\t\. 

There are two (and possibly many more) ways one can incorporate this choice within a calculation: 
1. Define the time integration in the time evolution operator to run over a contour in the complex plane: 
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With this choice Eqn. ((T|) becomes 



(5) 



Once this is done, complex conjugating the time evolution operator means that the time-forward contour does 
not coincide with the time-backward contour. The vacuum specification has broken the time-symmetry of the 
forward and backward time integrals. 

2. Analytically continue the interaction Hamiltonian occurring in the time evolution operator so that it becomes 
a function of a complex variable: 



H{t) ^ H{t{l + te)). 



(6) 



This is achieved by analytically continuing each of the times occurring in the expansion of Eqn. (JlJ . Since the 
momenta, fc, generically appear with the conformal time r in the combination kr this procedure is equivalent to 
an analytic continuation of the momenta flowing through the vertex. In practice one needs only to analytically 
continue the times or momenta appearing in the exponents of the mode functions. 



If we ignore this issue for a moment it is straightforward to express Eqn. ^ as [3, [37 



{w{t)) = 



■N 



N=0 



dt 



N 



dtN-1- 



dh {[Hintih), [i/i„t(t2), ... [HintitN), W {t)] 



(7) 



Eqn. (I7|) is formally consistent with Eqn. ([T]). However, this manipulation is only self-consistent if we use the second 
vacuum specification prescription above. Recall that for any symmetric, holomorphic function, /(ti,t2) — f{t2,ti), 



dU 



dt2f{tiM) 



dti 



dt2f{tl,t2). 



(8) 



Because / is symmetric under the exchange of its arguments, the integral over the square region on the left hand side 
of Eqn [8] is twice the right hand integral over the lower triangle, where t2 < ti. This result is easily generalized to 
more variables and moving from ([T]) to ([7|) requires repeated manipulations of this form. This manipulation requires 
that the integrations be interchangeable. The vacuum specification in Method 1 above breaks the equivalence of the 
integrals arising from operators on the right and left of W in Eqn. Consequently, terms like HWH, together 

with the contour specification, prevent one consistently writing down Eqn. ([7]), as the lower triangle is no longer 
identical to the upper triangle due to the asymmetry of the contour specification. Persisting with this approach risks 
losing information about part of the region of integration. In the second prescription, no contour is specified, so the 
manipulations above are perfectly safe. 

However, in applying Eqn. ^ one splits terms (e.g. H{ti)W{t)H{t2)) arising from Eqn. ([T]), which must be 
summed before any limits are taken, in order to avoid introducing unphysical divergences. At second order, Eqn. ([7|) 



(Wit)), 



dt2 



dti{{Hint{tl)Hint{t2)W{t) + W{t)Hint{t2)Hint{tl) 



-Hint{tl)W{t)Hint{t2) - Hint{t2)W{t)Hint{tl))) 



The terms appear to occur in conjugate pairs, so it may appear that this expression reduces to 



{w{t))2 = -2n 



dU 



dti {Hint{tl)Hi^t{t2)W{t) - Hint{tl)W{t)Hint{t2)) 



(9) 



(10) 



This manipulation is valid for HHW and WHH, but the vacuum prescription above prevents one from writing the 
HWH terms in this fashion. 

Finally, one might proceed from Eqn. ([7]) without employing either of the prescriptions above and regulate the 
oscillatory integrals in the far past by adding the appropriate small imaginary component to the initial point. In this 
case Eqn. (1101) is again self-consistent and for many diagrams this approach will work without difficulty. The in-in 
formalism associates each vertex with a time integration, and provided each vertex has a distinct momentum flowing 
through it, these distinct momenta effectively keep track of the distinct regions of integration. Adding all possible 
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permutations of the momenta which arise from the sum over contractions - as per Wick's theorem ^3§\ - sums over 
all tessellations of the restricted integration region and picks up all contributions. This approach fails when the sum 
of the momenta at each vertex are not distinguished from one another - as would be the case if one decides to make a 
specific choice about the external momenta before the time integrals are done - in which case it can lead to spurious 
divergences, as we will see below. Consequently, our preferred approach is to work directly with Eqn. ([5]) rather than 
Eqn ([7|). This not only avoids unphysical separation of the terms, but involves much less algebra. 



III. GRAVITATIONAL WAVES FROM LOOPS OF SCALARS 



We begin with the one-loop correction to the graviton (tensor) power spectrum generated by loops of second order 
scalar modes. In spatially flat gauge, the degrees of freedom are scalar field fluctuations 6(j) and transverse-traceless 
metric fluctuations (gravitons) jij . At leading order in slow roll S(j) couples to gravitons through the 3-point interaction 

m 

S=\j d^xdTa^{T)Y'dMdjS(l), (11) 
vanishes by gauge choice. This yields the interaction 



where t is the conformal time and the interaction Sijj 
Hamiltonian, which, after moving to Fourier space, is 



Hintir) - ifc'*fc"^7„(k,r)<5(/>(k',r)J0(k",r). 



Expanding the free fields in Fourier modes 



{2nf 



(2^ 



[a{k)U,{T)+aH~k)UUr)], 



^7.,,k = 7.,(x,r) = J2 y (2^e^'^-[6^(k)4.(k)7|(r)+6t^(-k)4.(-k)7r(r)] 



(12) 

(13) 
(14) 



S = +, X 



where the polarization tensors are normalized so that s' ^tj^^''''^ — 2^** and the mode functions, Ukir) and jk{T), 
are given by the solutions to the free field equations of motion (obtained from the quadratic part of the action, see 
for example [1^.) In the de Sitter limit these are: 



7.(r) = -^(l-zfcr)e*^-. 

The propagators are then 

{5Mr)5Mr')) = C/,(t)C/; (r')<53(k + p), 
(7,- ■,kW7/l,p(r')> - 7fe(r)7;(r')4(k)^L(-p)^.."^(k + p). 
To compute the one loop contribution to the tensor power spectrum, we evaluate^ 



(15) 
(16) 



(17) 
(18) 



dTi {Hint {Tl)Hint (t2 )7'^ {t* hi] (t*)) 



+ I dTi I dT2{Hint{Tl)Y\T^hi]{T*)Hint{T2)) ■ 

-OO-I- 



(19) 



^ The 1-vertex 1-loop scalar loop contribution to the graviton 2-point correlator is scale-free, and does not generate log corrections |23ll24(|. 



FIG. 1: The Feynman diagram corresponding to 2-point graviton correlation with a scalar loop. As we explain below, the 
diagram may be evaluated off-shell, so the momenta on the external legs have distinct labels, to avoid unphysical divergences. 



where —oo± = — oo(l ± ie) denotes the contour choice, the de Sitter limit a short calculation gives 

d^p' f d^p" 



7fc(T"*)7fc(T*) 



dry 



lk{T2)U;,{T2)U;„{T2) 



dri 



-\ll{T*)lk{T*) [ '^lk{T2)U*,{T2)U*„{T2) [ (t1 ) C/p' (ti ) C/p/- (n ) 

^ J -00+ '2 J-oo_ '1 



where 9 is the angle between the external momentum k and the momentum p'. However, if we use (jlOp to perform 
the calculation the third line above is replaced by (dropping irrelevant prefactors) 



5ft 



^l*k{r2)U;,{T2)U;„{T2) 

'^2 



— 7fe(n)J/p' (Tr)J7p"(Ti) 

^2 



-00 ' 1 



7fc(T*)7fc(r,) / —{l+iaiT2-a2T2-ia3T2)e' 
-00 '^2 



-iaiT2 



dr I 2 X N 

— 2-(l — iaiTi — a2Ti + ia^Ti )e 



3\^miTi 



00 '1 



(21) 



where ai — k + p' + p" , 02 = k{p" +p') +p'p" and 03 — kp'p". 

In this case if we do not carefully track the vacuum specification we will always run into a divergence. If one 
performs the integration on-shell, the lower limit in the first integral may be made to vanish by adding a small amount 
of evolution in imaginary time 10]. However, because the momentum flowing through each vertex is identical, the 
remaining integrand contains factors of e°'^'^ and e^"^^"^ which cancel, rendering the second time integral divergent. 
However, by writing down ()2ip we have implicitly ignored the difference in the contours between r and r'. With the 
vacuum specification explicitly included, the first integral leaves — i(ai(l + ie) — ai(l + ie'))T in the exponential in 
the outer integral.'^ Consequently as long as e, e' 7^ the integral is finite. The final limit e, e' — > is rendered finite 
and order-independent by symmetrizing over e, e', which is now required because H(ti)W H{t2) and H{t2)WH(ti) 
are no longer conjugates. Performing this symmetrization is equivalent to performing the integration in the opposite 
order and thus picks up the other term. Including the vacuum specification information in this fashion is equivalent 
to performing the calculation off shell and summing the result, before taking the on-shell limit. Pursing this approach 
one obtains a term of form 



7fe(T*)7fc(T*) 



2a2a3 2a3 



1 



ai 



(22) 



Much of this work can be avoided if one begins from Eqn. fS]), which yields Eqn. ((20|) . The integral factors into an 
expression of the form | J dr]f{'r])\'^ , which is not only well behaved everywhere but also reduces the double integral to 
a single integral. Carrying out the integration yields 



7fc(T*)7fc(T*) 2 



l(!ik^p''^p' 



//3 



Q3T 
ai 



i{aia2 + 03) 



037- 



(23) 



We are grateful to Xingang Chen |40|| . David Seery, Martin Sloth and Filippo Vcrnizzi for discussions about this point. 



FIG. 2; A sample diagram with N = 8 legs, two m — 8 point vertices and p = 4 internal lines. 



which is quickly shown to be identical to the expression in Eqn. (j22p . The full expression is 



(27^)' o n.r. f f d'p" P' . 4 



-ePcc(fc)^77(fc) 



(27r)3 J (27r)3 p"^ 



sin'' 6l5(k + p' + p") 



4fc3 + 2p'^k + 5^^ _ (6fc^ + 5/^) p" -6fc^ - bp'^k^ - Ap'^k - 5p'^ 
k k ^ kK 

3 (p"* + 2fcp'3 + + 2fcV) 



p"k^+p"^k 

X2 



(24) 



where K = k + p' + p" . The second and third lines come from the 1st term of Eqn. ((20| while the remaining terms 
come from the second term of Eqn. ((20|) . This expression takes the form 



(27r) 



eP^^{k)P^^{k) J d?p' sin^e j (fp"Y^ 



fa{p',p",k) 



(25) 



where fa denote functions of {p' ,p" , k) associated with the a-th power of K. We can eliminate the sin* 9 term in the 
integral with the useful identity 



j (fp j dV'5^(P + p' + k) sm^"ef{p,p'k) ^ 

'K{2k -K){K~ 2p'){K - 2/')"' " 



(-1)" 



27r 



pdp 



p+k 



p'dp' 



\p-k\ 



fip,p\k). 



(26) 



We are only interested in the logfc contributions to this integral, and the sin* 6 term effectively contains a factor K^, 
so the only terms in equation (j25[) that are not simply polynomial divergences are those with a > 2. After some 
straightforward algebra, we obtain 



dp' 



p'+k 



dp" 



p' — k\ 



pUp/u ^p'2p"k + Bp'-^p" + 5p"^p'^ + p'p' 



kK 



,(27) 



where the ellipses indicate polynomial terms which we will drop. Eqn. ([7f|) is IR finite, as the integrand contains no 
negative powers of p' or p" . Following p^ . we dimensionally regularize this expression to obtain 



35 



t^^PqC log k 



(28) 



where = H /{AM^k e) and P^^ — H /{M^k ), which are the uncorrected primordial curvature and tensor power 
spectra respectively. The loop corrections to the tensor power spectrum via a scalar loop are ^ (H/Mp)'^, and thus 
minute, as we would expect. 



IV. THE N-POINT 2- VERTEX INTERACTION 



The splitting of the terms in the expansion of ([T]) can induce spurious divergences in general multi- vertex topologies, 
as well as the 2-point diagram considered above. These divergences necessarily cancel in any careful calculation, but 
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we now show that the problem can arise for a general iV-point function. We consider the class of diagrams with N 
external legs, 2 vertices, and p internal lines - an example is given in Figure [3 

Suppose we have a term in the interaction Hamiltonian of the form /(r)^0™(T), where /(r) is some coupling. At 
second order one generates a general A'^-point function through a 2- vertex process with p internal lines. For simplicity 
we restrict ourselves to symmetric diagrams where each vertex has the same form, and with no loops which begin and 
end on the same vertex, so there are p = m — N/2 internal lines and N is even. At second order by this diagram is 
generated by 



= -23? 



J —OG J —OC 



dri/(ri)(5Ct(n)<50(n)~<5Ct(^2)) 



dT2f{T2 

If one had used the expansion of Eqn. ([7|), one would have 

dT2f{T2) r dri/(Tl)(<5Ct(Tl)<5Ct(^2)<50(T,)^) 



(29) 



dT2f{T2) 



dnfin) (<5C\(n)<50(r.)^<5C\(r2) + <5C\(^2)<50(r.)^5Ct(n)) . (30) 



where the second term in Eqn. has been split into two pieces in Eqn. ([50)1 . In the de Sitter limit the second term 
becomes 



N/2 



tt2 " tt2 

j=N/2 1 



P 

n 



d^k 

(27r)3 2ql 



N/2 



N/2 



(1 - z<7fcri)(l + *<7fcT2)e*«^(^^-^^)<53(^ k, - ^ qfe)'5'(E + Y1 "^^^ 



k=l 



k=l 



+{N\ — 1) permutations . 



(31) 



The product of delta functions ensures overall momentum conservation, the p\ counts the number of equivalent ways 
of contracting the internal lines, and the final sum over the permutations counts contractions of the external lines into 
the vertices. At this stage, the only difference between Eqn. (l29l) and Eqn. pO|) is the limits of integration, and that 
Eqn. (|30p requires that we add the same term with the order of integration reversed. The first term can be obtained 
from this expression simply by flipping the signs of the ki where i £ {1, A^/2}. Dropping all the irrelevant factors, 
suppressing the momentum integral and putting in the integration over the times. 



f{T2)dT2R{T2). 



N/2 



N 



/Ta 
/(ri)driQ(ri)e'(5:fii fc.+E. 
'CXJ 



'-5(^ - ^ qfc)5( ^ ki + ^ Qfe) + (A^! - 1) permutations. 



(32) 



k=l 



^N/2 

rN/2, 



k=l 



where Q and i? are the polynomials (3(ri) = ]\i=t{'^-iktTi)Y\!'i,{l-iqkTi) andi?(T2) = njljv/2(l+*^j'^2) Hf (l+^9i'^2)- 
The upper limit of the T2 integral is labelled by Tq, and we will take it to be either Ta — or Ta = T2. If we begin with 
Eqn. ([2^ . Ta = T, and Eqn. (|3ip is simply the product of two integrals. Moreover, after exchanging the momenta 
these integrals are simply complex conjugates of one other. To evaluate the term one must perform a single integral. 

Conversely, if we use Eqn pop , Ta — T2 and we encounter similar complications to those seen with the loop in the 
previous section. Performing the first integral in Eqn. (|32[) with Ta = T2 leaves 



/(T2)dT2i?(T2)5(T2)e^(^-i >'^~^"=-/2 '=-)-M(^ - ^ qfc)<5( ^ + J2 'f^'^ + " l)permutations, (33) 

1=1 fc=l j=N/2 k=l 



N/2 



N 
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where S{t2) is another polynomial. Now, if one restricts to X^i^i ~^ J2'j=N/2 % before^ performing these integrals 
one will encounter divergences. These divergences are spurious as discussed above. Properly including the vacuum 
specification prevents the limit from vanishing. Once the vacuum information is included the limit is really 'Y^^Si 
e'). Adding all the terms together (including those with the integration order reversed, which we have not written 
down) renders the result finite and order independent in the limit e, e' — > 0. 



Overall momentum conservation, enforces X^iL 



but does not require X^iLi ~ 



12j=N/2 % • Consequently, problems only arise if we render the external momenta indistinguishable before performing 
the integrals. For the special case where N = 2 overall momentum conservation already requires that the external 
momenta are identical, so this problem shows itself immediately when computing simple loop corrections to the 
propagator. Employing an explicit vacuum specification at the outset amounts to working off-shell by an infinitesimal 
amount, and thus ensures that the calculation is divergence- free. This problem can arise in any multivertex diagram 
and is not restricted to loop graphs - to see this, simply set p — 1 in our general topology, reducing it to an TV-point 
2-vertex tree diagram.^ 



V. CONCLUSIONS 



We have explored the role of initial state selection in the operator formulation of the "in-in" approach to quantum 
field theory, demonstrating that this requires careful handling at second order and above. Issues arise from the two 
different time contours which select the vacuum |in) and (in| at early times, and thus similar issues are not encountered 
in conventional in-out computations. To perform calculations, one uses the time evolution operator, U{t, tq) to evolve 
the "in" state forward in time to t, at which the correlation is to be computed. We show that many complications are 
avoided by keeping the initial state prescription explicit throughout the computation. Moreover, avoiding expansions 
which unnecessarily split terms into pieces which are individually unphysical simplifies many calculations and sidesteps 
these vacuum issues entirely. 

This issue was first encountered in 1-loop corrections to the scalar propagator [2^, and we again consider a 1- 
loop calculation here. We calculate the correction to the primordial gravitational wave spectrum produced by scalar 
bremsstrahlung during slow roll inflation. As one would expect, this is tiny, of order ~ eP-y-y(fc)P|^(;(fc) In(fc) where 
P^^(fc) and Pry^{k) are the (tree level) spectra of primordial curvature perturbations and gravitational waves respec- 
tively, and e is the usual slow roll factor. We show that one can avoid (potential) spurious divergences and messy 
algebra by working only with the physical terms arising from Eqn. ([T]) rather than the expansion in Eqn. ([7]) which 
subdivides some terms into pieces which, considered alone, are unphysical. 

We then show that one can encounter similar problems with these terms in any diagram with two or more vertices, 
including tree-level expressions. We demonstrate this by considering the A^-point function generated by a diagram 
with 2-vertices and p > I internal lines. Here, spurious divergences are not ubiquitous, but appear if one specializes to 
highly symmetric momentum configurations before performing the time integrations. This problem is avoided if the 
vacuum selection is explicitly and consistently imposed throughout the calculation, even in these specialized limits. 
In addition, eschewing seemingly convenient splits simplifies the resulting algebra, as a double integral is replaced 
with the product of a single integral and its complex conjugate. Finally, while our analysis involves only 2-vertex 
interactions, Eqn. ([7]) splits up terms at every order, so this issue can arise in any multivertex diagram. None of these 
divergences will be physical, and in all cases the problem can be ameliorated by imposing the contour choice at the 
beginning of the calculation. 

While it may seem that many higher order diagrams in cosmological perturbation theory are of purely academic 
interest as they relate to intrinsically miniscule effects (for an example where higher order effects can be important, 
see [4l'|), as the quality of astrophysical data and the theoretical sophistication of very early universe models increases, 
non-leading contributions to cosmological perturbations will become increasingly important. 



^ Faced with a complicated diagram, it is always tempting to restrict attention to situations with a high degree of symmetry. An obvious 
example is the equilateral limit where all the momenta are of equal magnitude, which is a common choice in calculations of the 3-point 
function. 

® With N = A and p = 1 the topology matches the graviton [3l| and |32| scalar exchange contribution to the scalar 4-point (trispectrum). 
While the specific details of the interaction differ, the algebra is very similar to that of the diagram we consider here. 
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